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I summarize recent progress in lattice gauge theory, with particular emphasis on results from numerical simulations. A major 
success has been the determination of the light hadron spectrum in the quenched approximation with sufficient accuracy to expose 
statistically significant disagreements with the experimental spectrum. The differences are, however, fairly small, ~ 5 — 10%. The 
data are also accurate enough to show evidence for artifacts of quenching predicted by chiral perturbation theory, f give an update 
on results for light quark masses, the kaon B-parameter, and the decay constants and B-parameters of heavy-light mesons. Most 
of these are known in the quenched approximation to ~ 10% accuracy or better, and preliminary estimates of quenching errors 
are of comparable size. One exception is the light quark masses, for which the quenching errors appear to be larger. I discuss the 
computational requirements for simulations of QCD with all approximations controlled, and argue that they will likely begin once 
computers sustain about 10 Teraflops. This is 30-40 times faster than present state-of-the-art machines. This estimate assumes that 
improvements in the discretization of lattice fermions are sufficient to allow continuum extrapolations to be made with a minimum 
lattice spacing of 0.1 fm. 1 review results obtained with improved discretizations and conclude that they satisfy this requirement in 
most cases. Examples of successful improvement include the calculation of the glucball spectrum and excited heavy-quark potentials 
in pure Yang-Mills theory. Finally, I discuss recent developments which may allow simulations of QCD with full chiral symmetry 
even at finite lattice spacing. 



1 Introduction 

In this talk I aim to give an overview of recent progress 
in lattice gauge theory, to report on the status of cal- 
culations of various quantities of phenomenological in- 
terest, and to give my expectations for future progress. 
I mainly consider results from numerical simulations of 
QCD, although towards the end I discuss some "analyt- 
ical" advances. Given constraints of time and space, I 
have been highly selective, and I apologize in advance 
to those whose work I have not included. An excellent 
source for more details is the proceedings of the annual 
lattice conference, the most recent of which was held in 
July 1998EJ. Previous reviews in the LCHEP and Lepton- 
Photan series (by Luscher last yearH, and by Flynn in 
19960) have different emphases than mine and are also 
highly recommended. 

The most important feature of the lattice is that it 
provides a non-perturbative regularization of field theory. 
This is necessary in order to investigate the properties of 
strongly coupled theories. The prime example is QCD at 
low energies or temperatures, and this is the major fo- 
cus of numerical lattice simulations. There is, however, 
another non-perturbative sector of the standard model, 
namely the electroweak sector close to the finite tem- 
perature symmetry-restoration transition. The proper- 
ties of this transition are being actively investigated us- 
ing lattice methods. These simulations are simpler than 
those for QCD, since fermions can, to a good approx- 
imation, be integrated out analytically. The status of 
calculations is described in the recent review of Laine 
and Rummukaincn and I will not discuss them here. 

The ultimate goals of simulations of QCD include 



the following: 

• To test that QCD is indeed the correct theory of 
the strong interactions; 

• To make predictions for exotic mesons and baryons; 

• To determine the host of hadronic matrix elements 
which are needed to extract the underlying quark 
parameters (masses and weak mixing angles) from 
experimental results. 

These are the topics I discuss in this talk. There are 
other important goals which I will not discuss: 

• The attempt to gain at least a qualitative under- 
standing of the physics of confinement and chi- 
ral symmetry breaking (a recent review is that by 
NegeleH); 

• To understand the behavior of QCD at finite-tem- 
perature and density (see reviews by KanayaB and 
LaermannQ). 

Another interesting topic I will cover only in passing is 
the status of attempts.iQ_discretize chiral gauge theories. 
For reviews see Refs.Ejtffil 

The outline of this review is as follows. Following a 
lightning description of lattice calculations, which serves 
also to introduce some common terminology, I explain 
why simulations of lattice QCD are so computationally 
intensive, and indicate how close we are to controlling 
all sources of error. In sec. |^ I describe the status of 
quenched calculations of the spectrum. Section ^ gives an 
update on the successes of the "improvement program" 
at reducing discretization errors, successes that are illus- 
trated using the examples of the spectrum of glueballs 
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and potentials in the SU (3) Yang- Mills theory. Section [| 
gives a brief update on the status of spectrum calcula- 
tions including dynamical quarks. I then give a summary 
of the results for those hadronic matrix elements which 
are calculated most reliably. Section [?] gives an overview 
of the recent progress in formulating lattice QCD with 
exact chiral symmetry. 

2 Why are lattice QCD calculations hard? 

Numerical simulations of gauge theories have been un- 
derway since Creutzls seminal work almost twenty years 
ago (see, e.g., Ref. Ell). The first simulations including 
fermions (in. the quenched approximation) followed soon 
thereafterB. Yet we are still not at the point of having 
reliable results for QCD. In this section I will describe 
why the simulations are so difficult, and explain what it 
will take to get reliable results. 

I begin with an quick review of lattice methodology. 
Wilson taught us how_to discretize QCD while maintain- 
ing gauge invariance E3: fermion fields are placed on the 
sites of a lattice, while the gauge fields live on the inter- 
vening links. The simplest lattice, used in most numerical 
simulations, is a hypercubic array, with equal spacings in 
all directions. The symmetry between sites makes the 
system well suited for numerical implementation on par- 
allel computers. It is advantageous for some applications 
to make one or more lattice spacings different, and I will 
mention an example of this below. For numerical simula- 
tions one must also truncate to finite volume. What re- 
sults is a finite dimensional quantum mechanical system 
that can be simulated using Feynman's sum over paths. 
This is only practical if the time direction is imaginary, 
so that one is in effect studying 4-dimcnsional statistical 
mechanicals. 

The fermionic path integral is over anticommuting 
variables, and is usually done by hanoQ, resulting in the 
partition function 

Zqcd = f[dU}cM-S g )Y[det{put + m^ n ). (1) 
J g 

In other words, the path integral is over all gauge fields 
U (one SU(3) matrix per link), weighted by the gauge 
action S g and the effective action resulting from fermion 
loops. S g is a discretization of the continuum gauge ac- 
tion—and is commonly taken to be that proposed by Wil- 
son E3. This has discretization errors of 0(a 2 ). pi at is 
a discretization of the Dirac operator. The most com- 
mon choices are Wilson fermions £3 (or improvements 
thereof which have_smaller discretization errors), and 
staggered fermions c3. Unimproved Wilson fermions have 

a There is, howevfpqa recent paper in which the fermion integral is 
done numericallylifl. 



discretization errors of 0(a), whereas staggered fermions 
have errors of 0(a 2 ). 

Physical quantities are obtained from correlation 
functions. The simplest examples are two-point func- 
tions such as that of the pseudoscalar density, 

C P {x) = -(u-ysd(a;) d 75 u{0)) . (2) 

From this one can extract the mass of the ir~ (and, less 
easily, of its excited states). Integrating out the fermions 
analytically gives the following expression 

C P (x) = Z Q i D f [dU] exp(-5 fl ) J] det(D lat + mf ") 

xTr r(^^)./ 5 (^T^)„J ' (3) 

which is suitable for numerical evaluation. Thus, to cal- 
culate a correlation function, one has to calculate quark 
"propagators" (i.e. inverses of the discretized Dirac op- 
erator plus mass term) in the background gluon config- 
uration, connect them together as shown by the trace, 
and average the result over gauge configurations weighted 
by the total effective action. The average is done using 
Monte-Carlo methods. 

A very important distinction in the following is be- 
tween "dynamical" and "valence" quark masses. The 
former are those which appear in the determinant, while 
that latter occur in the propagators. I have labeled the 
quark masses in Eq. |j| according to this terminology. In 
QCD, these two masses are the same, flavor by flavor, 
since they both arise from the same original action. But 
once one has integrated out the fermions, one has the 
option of picking the masses to be different. This is 
called the "partially quenched" approximation, and is 
made in most large scale calculations at present. Typi- 
cally, the dynamical quark masses are heavier than the 
valence quark masses. The fully quenched limit is that 
in which m dyn — > oo, in which case the fermion deter- 
minant is a constant and drops out of the action. In 
summary, the partially quenched approximation means 
taking m val ^ m dyn , and the limit m dyn — > oo with finite 
m val is the quenched approximation. 

It is important to realize that (partially) quenched 
theories are unphysical — for example, they do not have 
an hermitian Hamiltonian. These approximations are 
made simply out of necessity because of the difficulty 
in simulating QCD with light dynamical masses. To il- 
lustrate this difficulty I will use an approximate scaling 
law for the number of floating-point operations (Flops) 
that are required to generate a statistically independent 
gauge configuration using 2 flavors of_dynamical stag- 
gered fermions and present algorithms E3 

# Flops _ 6000 x N^N t 
Indep. Config. ~ (aTO dyn ) 5 / 2 
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Table 1: Time for 100 QCD configurations on a Terafiop machine. 



m dyn 


ra n L 

m P fm 


a=.2fm a=.lfm a=.05fm 


m s /2 


0.56 2.5 


7.5/mins 11/hrs 42/ days 


m s /8 


0.3 3.5 


15/hrs 58/ days 1.4/ yrs 


m s /25 


0.18 6.7 


0.4/ yrs 3.7/ yrs 330/ yrs 



Here N s and N t are the number of points in the spatial 
and Euclidean time directions. Assuming a typical value 
for their ratio, N t /N s = 2, this converts to 



# Tflop-hours 
100 Configs 



11/ 



L 



2.5fm 



0.1fm\ 
a J 



6.5 



m s /2 



2.5 



(5) 

where L — N s a is the spatial size of the lattice in phys- 
ical units, and m s is the physical strange quark mass. I 
choose 100 independent configurations in the denomina- 
tor since this is roughly the minimal number needed to 
obtain small statistical errors in masses. I have intro- 
duced a factor / to account for differences in simulation 
times between different discretizations of fermions. Nu- 
merical examples of this formula are give in Table 0. Note 
that the lattice size is chosen to satisfy both L > 2.5 fm 
and m^L > 4.5 to avoid significant finite volume effects. 

The top row in the Table corresponds approximately 
to the lightest quark masses used in present simulations 
of (unquenched) QCD. Simulations at this quark mass 
are being done with the lattice spacings down to 0.1 fm. 
The bottom row corresponds to the physical average 
light-quark mass. This is not quite QCD, since a dy- 
namical strange quark is not included, but it gives an in- 
dication of the CPU time required for a direct simulation 
of QCD. Finally, the middle row corresponds to a quark 
mass which should be light enough that one can extrap- 
olate from it to the physical values. In other words, we 
probably do not need to use masses smaller than those of 
the middle row, because we can use chiral perturbation 
theory to guide us in doing extrapolations. The chiral 
expansion parameter is m^/m^, and this is 10% for the 
middle row (but 30% for the top row). 

A striking feature of these numbers is how it be- 
comes rapidly more difficult to simulate as the lattice 
spacing and quark mass decrease. This follows directly 
from Eq. [| Part of this increase is unavoidable (e.g. 
the number of lattice points increases as a decreases) but 
part is due to shortcomings of present algorithms (e.g. 
the critical slowing down as m q — > 0). Although I expect 
algorithms with improved scaling laws to be developed, 
I use the scaling of present algorithms for my estimates. 

The dependence on a makes clear the importance of 



the "improvement program" , which aims to reduce dis- 
cretization errors and thus allow simulations at larger lat- 
tice spacings. In quenched simulations with unimproved 
Wilson fermions, it turns out (as discussed below) that 
one needs to go down to o k 0.05 fm in order to extrapo- 
late to the continuum limit. It would be a great advance 
to be able to work with a minimum lattice spacing of 
a « 0.1 fm instead. I will show later that this is likely 
to be possible using improved Wilson fermions. With 
staggered fermions, on the other hand, it appears that 
a w 0.05 fm will be necessary. 

Let us assume, then, that we work with improved 
Wilson fermions, and can get away with a minimum a of 
0.1 fm, and a minimum quark mass of m s /8. Then we 
need only consider the top left 2x2 block of the table. 
What does this imply for the absolute time required? 
To answer this I must discuss the factor /. This repre- 
sents the extra cost of doing simulations with improved 
Wilson fermions compared to staggered fermions. .-Based 
on the experience of the CP-PACS collaboration 113 this 
factor is at least 10, implying that the simulation with 
m dyn = m s /8 and a = 0.1 fm would take almost 2 years 
on a sustained Terafiop machine. To be realistic one 
should also include a sizeable additional factor to account 
for the extra cost of including the dynamical strange 
quark, the need to do multiple simulations in order to 
extrapolate, etc. Thus I expect that reliable results for 
QCD itself will begin to emerge once we have machines 
sustaining 10 Teraflops.^] I want to stress that this is 
only 30-40 times faster than the present state-of-the-art 
machines. These are the CP-PACS machine at the Uni- 
versity of TsukubaJapan, which peaks at 0.6 Tflops and 
sustains 0_3TflopsEj; and the two Columbia-Riken/BNL 
machinesEJ which together have a peak speed of 1 Tflop, 
and present efficiencies of about 25%. CP-PACS has 
been in production mode for nearly 2 years, while the 
Columbia machine (0.4 Tflop peak) has been in produc- 
tion mode for 6 months, and the Rikcn/BNL machine 
(0.6 Tflop peak) is just becoming fully operational. 

Figure [j] gives a slightly different view of the status 
of present simulations. It shows a two-dimensional sec- 
tion through parameter space (and thus hides the depen- 
dence on the lattice spacing and the strange quark mass). 
Theories along the diagonal are physical, while those off 
the diagonal are partially quenched and thus unphysical. 
The fully quenched theory is on the righthand vertical 
axis at m dyn — oo. The upper shaded region roughly 
indicates where present simulations are done. Note that 
it is relatively inexpensive, for a given dynamical quark 
mass, to map out a range of valence quark masses. Thus 
progress will occur along the various arrows pointing to- 



b This estimate is somewhat more pessimistic than that I gave in 
the talk. Simulations with improved actions are somewhat slower 
than I had assumed. 
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Figure 1: Approaching QCD. 



wards the physical point in the LH corner. I wish to em- 
phasize that, because of the ease of collecting partially 
quenched data, it would be very useful if we could use the 
information from these unphysical theories to aid the ex- 
trapolation to the physical theory. This can perhaps be 
done using chiral perturbation theory, and work in this 
direction has begunElEi There is, however, a limitation 
on how small one can take m val , because certain artifacts 
of quenching diverge as m val — ► 0. Examples of this in 
the quenched theory will be shown below. 

To conclude this section, I list the important issues 
that effect the progress towards simulations of full QCD. 



• How large are quenching errors? To what extent 
can we use quenched QCD as a guide to QCD? 

• How small a lattice spacing is needed to extrapolate 
to the continuum? By how much can we increase 
a m i n with an improved action? 

• How small a quark mass is needed to reliably ex- 
trapolate to QCD? 

• Can we use (partially quenched) chiral perturba- 
tion theory to aid extrapolation to the physical pa- 
rameters? 

3 The spectrum of quenched QCD 

The CP-PACS collaboration presented this year final re- 
sults from their study of the spectrum of QCD in the 



Tabic 2: Parameters of CP-PACS quenched simulations. 



P = 6/ff 2 


a (fm) 


size 


La (fm) 


Configs. 


5.90 


0.10 


32 3 x 56 


3.3 


800 


6.10 


0.078 


40 3 x 70 


3.1 


600 


6.25 


0.064 


48 3 x 84 


3.1 


420 


6.47 


0.05 


64 3 x 112 


3.2 


150 



quenched approximation ESEfl. Their aim is to control all 
errors, except that due to quenching, and thus to obtain 
definitive results for the quenched spectrum. They use: 

• the simplest discretizations of the gauge and 
fermion actions, both due to Wilson. Thus the 
dominant errors are of O(a), coming from the 
fermions; 

• a fixed physical volume large enough for finite size 
effects to be negligible for single particle masses; 

• four (small) lattice spacings to extrapolate to the 
continuum limit (a = 0); 

• five valence quark masses, spanning the range 
m^/rap — 0.4 — 0.75, to extrapolate to the physical 
light quark masses; 

• large ensembles of gauge configurations so that the 
final errors, after extrapolation, are of size 1 — 3%. 

Some parameters of their simulations are collected in Ta- 
ble [| I focus mainly on the CP-PACS results for the 
spectrum since their computational power outstrips that 
of other groups by an order of magnitude. 

This calculation supersedes that^ione using the 
GF11 machine and presented in 1994 EEl. The improve- 
ment, due to more than an order of magnitude increase 
in computer power, is evident from the final results, 
shown in Figs. [2] and ||. These figures give the masses 
of pseudoscalar and vector mesons and octet and de- 
cuplet baryons. To understand them it is useful to re- 
view how the parameters of lattice simulations are fixed 
in quenched simulations. First, for each lattice spacing 
one extrapolates in the average light quark mass to the 
value at which the physical m n /m p ratio is obtained^ 
The resulting rho mass is then used to determine the 
lattice spacing (one extracts m p a from the lattice, so 
a = (m p a)/mjj hys ). At this point all masses can be con- 
verted to units of GeV (although one should keep in mind 
that one is really determining the ratio of the masses to 
m p ). Next, one interpolates in the lattice strange quark 

c Note that isospin breaking effects are being ignored. 
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Figure 2: Quenched spectrum with m s fixed using nig. 



Figure 3: Quenched spectrum with m s fixed using 



mass until a particular strange particle has its physical 
mass — in the present cases the K and c/> are used. And, 
finally, one extrapolates in the lattice spacing to a = 0. 
The results of this procedure are shown in the Figures — 
the pion and rho masses are not shown since they are 
input parameters. 

I draw two important lessons from these results. 
First and foremost, the quenched spectrum disagrees 
with experiment! This is an important check on the cal- 
culations, since we know that we are not calculating in 
the correct theory. The most significant disagreement is 
for the hyperfine splitting, — m^-, which is ~ 10% 

too small for either choice of m s . This is a 5 — 6a dis- 
crepancy. The nucleon mass is 7% too low (indepen- 
dent of m s ), although this is less significant (2.5a). Fi- 
nally, the splittings within the octet and decuplet are 
too small, although this discrepancy is greatly reduced if 
m s is fixed using m$. However, the Gell-Mann-Okubo 
formula for octet baryons and the equal-spacing rule for 
decuplet baryons are both well satisfied for either value 
of m s . 

The second lesson is that quenched QCD yields a 
spectrum which is reasonably close (5 — 10%) to that 
expected from QCD itself. In other words, although 
quenching is a drastic approximation, it gives us a model 
in which the gross properties of light hadrons are numeri- 



cally close to those in QCD. This is important because we 
want to have an idea of the accuracy of quenched results 
for other quantities, such as the hadronic matrix elements 
discussed below. One should keep in mind, however, that 
there is no universal size to quenching errors. 

Extrapolations play a key role in the CP-PACS anal- 
ysis, and I want to spend a little time discussing their reli- 
ability. I focus on the chiral extrapolation here, deferring 
discussion of the continuum extrapolation until the next 
section. To do the chiral extrapolation requires theoret- 
ical input as to the functional form, and this is provided 
by chiral perturbation theory (xPT). For the quenched 
theory one must use the quenched variant, Q%PT, and 
this yields the following prediction jfcjr_| hc chiral expan- 



sions of pion and nucleon masses E 

m, = 2Am 
M 



m N 



q (1 - S\og(2Am q /A 2 x ) + 0(m q )) 
SC imn + C 2 ml + C 3 ml + 0{m%) 



(G) 
(7) 



In these expressions, 6, A, A x , Mq and the Cj are con- 
stants which are related to the coefficients appearing in 
the quenched chiral Lagrangian. The result for m p has 
the same form as that for mjv but with different con- 
stants c3. The important features of these expressions 
are as follows. 

• Terms involving the constant S are artifacts of 



5 



quenching — they are not present in the chiral ex- 
pansions in QCD. They arise from loops containing 
the rf , which, in the quenched theory, remains de- 
generate with the pion. They give the dominant 
corrections to the leading order term as m q — > 0. 
In particular, the 8\ogm q term in ([3]) diverges in 
the chiral limit, although the contribution to 
itself vanishes in this limit. 

• Terms not involving 6 are present in both quenched 
and unqucnchcd theories, although the detailed ex- 
pressions for the coefficients differ. 

• Particle masses are predicted to be non-analytic 
functions of the quark mass — the non-analyticity is 
proportional to m q log m q in (||) and to m, ~ m X q 2 



and 



3/2 



In the past, most analyses have extrapolated using low 
order polynomials in m q , i.e. ignoring possible non- 
analytic terms. These parameterizations have adequately 
represented the data. The small statistical errors and 
long chiral level-arm of the CP-PACS calculation make 
possible a test of the more complicated forms predicted 

QxPT. 

The most dramatic effect predicted by QxPT is for 
the pion mass: the quantity y = / (mi + nyA should 
diverge logarithmically at small m q . (See Ref. cB for the 
full expression for non-degenerate quarks.) This is to be 
compared to QCD, where y linearly approaches a con- 
stant in the chiral limit. Evidence for this effect has 
previously been_seen with staggered fermions (and is re- 
viewed in Ref. but not with Wilson fermions. The 
CP-PACS results are shown in Fig. ||. The predicted rise 
at low quark masses is seen for all four lattice spacings, 
with the fits yielding consistent values of 8 ~ O.IfL±0.02. 
This rise has also been seen by the FNAL group Ell, who 
are able to work with very small quark masses by using a 
modification of the quenched approximation, and-by the 
QCDSF group using improved Wilson fermions cHl. The 
values for S are all roughly consistent with each other, 
and with those obtained using other quantities EJ. 

The most economical explanation of these results is 
that they are quenched artifacts due to 7/ loops, as pre- 
dicted by QxPT. This interpretation is supported by the 
fact that fits to the predicted chiral behavior of vector 
mesons, baryons and pseudoscalar decay constants (in- 
cluding those for particles composed of nprjjilegenerate 
quarks) represent the CP-PACS data wello'Ell. It is re- 
assuring that the lattice results are good enough to pick 
out physical effects coming from loops of light particles, 
and thus to make contact with the analytical predictions 
of QxPT. It also gives credibility to attempts to use 
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Figure 4: CP-PACS results for y at all four lattice spacings. Points 
with open triangles and squares are for non-degenerate quarks; oth- 
ers are for degenerate quarks. The quark masses are in the MS 
scheme at fi = 2 GeV. The lines are drawn to guide the eye and are 
not fits. 



QxPT to estijjriaJte the size of quenching errors in other 
quantities E^&Ea. 

On the other hand, one cannot rule out that the 
rise in y is due to a collusion between terms analytic in 
m q of linear, quadratic and cubic order. CP-PACS have 
investigated this possibility by fitting all the masses to 
polynomials of the lowest order required to give a good 
fit. The resulting light hadron masses differ at finite lat- 
tice spacing from those obtained using the QxPT forms, 
but it turns out that they agree within errors after ex- 
trapolation to the continuum limit (see Fig. [5] below). In 
particular, the hyperfine splitting is still significantly be- 
low the experimental value. Thus the fact that quenched 
QCD disagrees with experiment appears to be robust. 

The importance of understanding the chiral extrap- 
olations has also been stressed by the MILC collabora- 
tion, whp_^have done extensive studies using staggered 
fermions E3. 

My overall conclusion concerning chiral fits from 
present data is that, while it is clear that the lowest or- 
der polynomial fits are inadequate once one goes below 
m^/trip < 0.5, a complete resolution of the issue probably 
requires smaller masses than those used by CP-PACS. 
This is part of the reason why I chose m^/mp = 0.3 
(rather than the minimal value m^/mp = 0.4 used by 
CP-PACS) as the number to aim at in Table [l]. Actually, 
sorting out the chiral extrapolation may be easier in the 
unquenched theory, since there are no "enhanced" loga- 
rithms. On the other handj-such logarithms are present 
in partially quenched QCDcilH (and give rise to the ar- 
tifacts in the lower shaded region in Fig. |l|). 
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While the CP-PACS results are very impressive, 
they do rely on certain assumptions about extrapola- 
tions. It is clearly important to check them using at 
least one other fermion discretization, for which the sys- 
tematics will differ. The MILC results with staggered 
fermions, while not as extensive (e.g. they do not con- 
sider baryons composed on non-degenerate quarks), do 
provide a partial check. They quote the extrapolated 
value ui n I Tfip — 1.254 ± 0.018 ± 0.028, consistent with 
the experimental value 1.22. This is higher than the CP- 
PACS result of 1.143 ±0.033, but only by about 2.5cr. It 
would be very interesting to see results from staggered 
fermions for the K — K* hyperfine splitting. There are 



also a number of large scaL 
improved Wilson fermions 



itions underway using 
3, but these do not yet 
have a range of lattice spacings or statistical accuracy 
comparable to that of the CP-PACS collaboration. 



4 The improvement program 

The discussion of sec. || highlighted the importance of re- 
ducing discretization errors. This is essential if we are to 
have unquenched simulations of QCD in the near future, 
and it has been the focus of much research in the last 
few years. In this section I report on the progress in this 
direction. 

With an unimproved action, such as that for Wilson 
fermions, the discretization errors start at O(a) 



M^[l + «A 1 + (aA 2 ) 2 + 



(8) 



The size of the corrections is set by non-perturbative 
physics, so we expect Ai ~ A 2 ~ Aqcd- Examples of 
continuum extrapolations from the CP-PACS data are 
shown in Fig. ||. The fits are to a linear form (A 2 = 0), 
and give values for Ai in the range — 300 MeV. Note 
that a linear fit is reasonable for such small lattice spac- 
ings. If one introduces a quadratic term of the expected 
size the extrapolated values change by less than la. 

It is useful to estimate the minimal lattice spacing 
(ttmin) needed to extrapolate to the continuum with a 
given accuracy. I take the accuracy required after ex- 
trapolation to be 1%. To approximate the extrapolation 
procedure I assume that Mh a d and dM^d/da can be de- 
termined accurately at a m ; n , and used to extrapolate lin- 
early to zero. For the quadratic function in Eq. || the 
resulting error is — (a m i n A2) 2 . Setting this to be 1%, and 
taking the fairly conservative value A 2 = 400 MeV, I find 
o-min = 0.05 fm. This estimate is consistent with what 
CP-PACS finds in practice. 

Now assume that the calculation has been improved 
so that Ai = 0. The functional form is then quadratic 
+ cubic. Proceeding in the same way, except this time 
extrapolating linearly in a 2 , the error is — (a m i n A3) 3 /2. 
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Figure 5: CP-PACS continuum extrapolations, with m B fixed us- 
ing rax, and the scale in GeV set by m p . Solid (open) points 
are from chiral extrapolations based on QxPT (polynomial) fits. 
Experimental values are shown by diamonds at a = 0. 



Setting this to be 1%, and taking A 3 = 400 MeV, I find 
flmin = 0.14 fm. To be conservative I assume that there 
are some quantities with larger discretization errors (and 
there are known examples where this happens, such as 
Bk with staggered fermions), and so reduce the minimal 
lattice spacing to 0.1 fm. Even so, this doubling of the 
minimal lattice spacing compared to the unimproved case 
makes a huge difference in our ability to simulate QCD, 
as discussed in sec. |[ 

How does one go about improving the results? Most 
of the work so far has followed the Symanzik programs 
of improvement order by order in a. Symanzik showed 
that, without improvement, the long-distance (p<l/a) 
modes of the lattice theory can be described by an ef- 
fective continuum Lagrangian containing all operators of 
dimension 5 allowed by the lattice symmetries. This is 
analogous to the more familiar effective theories such as 
the Fermi theory of the weak interactions. For QCD this 
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reads 

£ off (a) = £ QCD + c 5 {a)aqcj^F^q + 0{a 2 ) (9) 

where is the gluon field strength, and C5 only de- 
pends logarithmically on a. There are other dimension-5 
operators but these can be absorbed by a chajage of vari- 
ables, and so do not effect physical quantities!^, ft is the 
C5 term which leads to the O(a) corrections in hadron 
masses discussed above. To remove these corrections, 
one must add to the lattice action the corresponding op- 
erator ell 

£\Z P = £ ?at + csw(3 2 )(ag<V,iVg)i at ( 10 ) 

and adjust the coefficient csw until C5 vanishes.^] If this 
single tuning can be done, then there are no linear dis- 
cretization errors in any physical quantity. Note that 
there are no gluonic operators of dimension 5 — they first 
appear at dimension 6 and give rise to errors of 0(a 2 ). 
Thus for the first stage of improvement for QCD the 
gluon action does not need to be improved (although it 
may still be advantageous to do so). 

The key then, is to have a method for determining 
csw- This can be done in perturbation theory, order 
by order in g 2 , and has been carried out in practice up 
to one-loop order. This leads to what are called "per- 
turbatively improved Wilson fermions" ,^| which have dis- 
cretization errors of size a 2 a and a 2 . It is clearly prefer- 
able, however, to have a non-perturbative determination. 
A method for doing so has been developed and applied 
by the ALPHA collaboration in the last few years, both 
for the quenched and two-flavor unquenched theories. I 
will not discuss the method since it was reviewed last 
year by Martin LiischerB. What is new this year are the 
results of extensive numerical tests of "non-perturbative 
0(a) improvement" as it is called. In other words, with 
csw determined non-perturbatively, do hadron masses 
behave as in Eq. || with Ai = 0? The answer is yes! 

4-1 Numerical tests of improvement 

Studies involving several lattice spacings have been done 
by the SCRI group BN, the QCE1SF collaboration E3, 
and the APFiTOV collaboration B, and the ALPHA 
collaboration a. I will show one plot, Fig. ^, which 
nicely captures the main features of all the results c3. 
To test improvement one need not extrapolate to physi- 
cal light quark masses — in fact it is better to use quarks 

d This is not as simple as cgw = — C5 because the addition of the 
extra term leads to a shift in C5 due to linear and higher divergences 
in lattice integrals. 

e Other names for the improved fermion action are the 
"Sheikholeslami-Wohlert" action ("SW" action for short) and the 
"clover" action (based on the form of the lattice version of F^y 
appearing in the improvement term). 
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Figure 6: Results for m p versus a 2 at fixed m^jmp = 0.7 in the 
quenched approximation. 

with roughly the strange quark mass which can be sim- 
ulated directly and for which the statistical errors are 
small. This is also advantageous because it avoids so- 
called "exceptional" configurations discussed below. The 
plot shows results for the vector meson mass my versus 
a 2 at a fixed value of the pseudoscalar to vector meson 
mass ratio. The lattice spacing is determined using the 
"string tension" a, which is the coefficient of the linear 
term in the potential between a heavy quark and anti- 
quark in the quenched theory. All we need to know about 
a is that it is a gluonic quantity and thus approaches 
the continuum limit with errors of Ola 2 ). Any linear 
dependence on a is due to discretization errors in the 
fermion action. The absolute scale is only approximately 
defined since this is a quenched calculation, but the ar- 
row in the figure labelled u j3 = 6.0" roughly corresponds 
to a = 0.1 fm. 

The results in the plot are from several actions in- 
volving different degrees of improvement. The most 
interesting comparison is between unimproved Wilson 
fermions (crosses and squares) and non-perturbatively 
improved Wilson fermions (circles). The former have 
errors of 0(a) while the latter should have 0(a 2 ) er- 
rors. The solid lines are the result of fits to these hy- 
potheses with the constraint that the continuum value 
be the same. The fit is reasonable, showing that non- 
perturbative improvement works. The quantitative suc- 
cess is also striking. For example, at a « 0.18 fm the 
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discretization errors in hadron masses are reduced from 
20 — 30% with Wilson fermions to 1 — 3% with non- 
perturbative improvementcj. 

If this data alone indicated the efficacy of improve- 
ment in general then one would be able to work at larger 
minimal lattice spacings than my conservative estimate 
of a m i n = 0.1 fm. The success is not uniform, however. 
HeitgerCj, using the non-perturbatively 0(a) improved 
action, finds that, at a — 0.1 fm, the discretization er- 
rors are only a few percent in most quantities, but that a 
quantity related to f p differs by 28% from its continuum 
value. 

The remaining points in the figure illustrate other 
interesting results. One can attempt to improve the 
gauge action so as to reduce the a 2 errors (which are 
large in certain glueball masses as we will see below). 
This improvement has only been implemented perturba- 
tively, reducing the errors to a s a 2 . One can then non- 
perturbatively improve the fermions on these new back- 
ground gauge fields. This should remove the 0(a) errors, 
but one should not expect much reduction in the a 2 er- 
rors compared to the usual gauge action, because one has 
not reduced the a 2 errors coming from the fermions. This 
is what is found, as shown by the diamonds in the figure, 
with the dashed line being the fit. In this case there is 
no advantage gained from improving the gluon action. 

The set of four points with the "burst" symbol (and 
no fit line) show the result of perturbatively improving 
the fermion action. Here one expects that there should 
be an 0{a) term, although reduced from that for unim- 
proved Wilson fermions by a factor of a s . Qualitatively 
this is what is found. A fit to just an a 2 dependence 
(a linear fit on the plot) does not extrapolate close to 
the correct continuum value. Note also that the er- 
rors at a ~ 0.2 fm are substantially larger than for non- 
perturbatively improved fermions. 

I conclude that, for almost all of the quantities stud- 
ied so far, Symanzik improvement not only removes the 
leading correction, but also reduces the size of the dis- 
cretization error substantially for the lattice spacings of 
interest. This is a very encouraging result. 

Staggered fermions provide an interesting exception 
to this rosy picture, as shown by the two uppermost 
crosses in the figure. Staggered fermions are already 
improved — the errors start at 0(a 2 ) — but it appears that 
the coefficient A 2 is very large. This is perhaps due to 
the fact that the effective lattice spacing is twice as large. 
Atte mpts tn improve staggered fermions are just begin- 
ning cj'EEtEj — and will be needed if they are to present a 
competitive alternative to improved Wilson fermions. 

What I have discussed so far is improvement of the 
S-matrix: masses and scattering amplitudes. We would 
also like to improve hadronic matrix elements, such as 
Jb- This requires improvement of operators as well as 



the action, e.g. for fg one needs to improve the axial 
current. This is not an easy task, since for each opera- 
tor there are new improvement constants, and one wants 
to have a non-perturbative method for determining all 
of these. So far, methods have been proposed for all 
bilinear operators, and have been implemented in sev- 
eral casesH'E 2 ]. Extension to more complicated operators 
should be possible. 

What is needed now is to extend the improved calcu- 
lations towards the chiral limit. In particular, one would 
like to do the improved analogue of the CP-PACS cal- 
culation, including chiral and continuum extrapolations, 
in order to check that the results agree with the CP- 
PACS benchmark, and to determine what is the mini- 
mal lattice spacing needed. In particular, can we use 
Omin = 0.1 or larger, as required to approach QCD in the 
near future? A number of calculations are underway (as 
reviewed by Kenwaya). An unfortunate obstacle is the 
presence of exceptional configurations. These are con- 
figurations on which the operator pi a t + m has a zero 
eigenvalue for a physical quark mass, and thus the prop- 
agator diverges at this mass. In the continuum these 
zero eigenvalues move to m q — and are the expected 
zero-modes associated with topologically non-trivial con- 
figurations. They are present at positive m q on the lat- 
tice because of explicit chiral symmetry breaking. For- 
mally, this makes the quenched approximatiopJll-defined 
at small quark masses with Wilson fermions ell. It turns 
out that this problem is exacerbated by non-perturbative 
improvement, and makes it difficult to study the effect of 
improvement on the quenched light hadron spectrum for 
lattice spacings larger than 0.1 fm. This obstacle may -he 
avoided using alternate forms of the improved actionem. 
It is important to realize, however, that exceptional con- 
figurations are quenched artifacts, and will not present 
a problem in unquenched QCD because they are sup- 
pressed by the fermion determinant. Thus they are a 
temporary obstacle. 

4-2 Other approaches to improvement 

I have concentrated so far on removing as completely as 
possible the corrections linear in a. Other approaches to 
improvement are also being pursued. Lepage has argued 
for pushing Symanzik improvement to higher order, us- 
ing tadpole-improved perturbationJheory to fix the ad- 
ditional improvement coefficients lj. Tadpole-improved 
perturbation theory has proven successful at predicting 
the dominant lattice corrections in many cases E3. The 
hope is to work on coarser lattices, perhaps as coarse as 
a = 0.25 - 0.4 fm, which implies 1/a = 0.5 - 0.8 GeV. 
For example, Alford et al. o argue that errors in hadron 
masses could be reduced to a few percent even at a = 
0.25 fm if one combines tadpole estimates of the higher 
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order improvement coefficients with a non-perturbative 
determination of the leading improvement coefficient 
cgw- Their numerical results, using the D234c action, 
for which the tree-level errors begin at a 4 , support this 
claim. 

I think that this approach is very interesting, but I 
am uncomfortable at the reliance on perturbation theory 
at such large lattice spacings. At a = 0.4 fm, for example, 
the relevant coupling constant is a s (Tr/a) w 0.4. This in- 
troduces two problems. First, one is essentially giving 
up on doing an extrapolation in a, since several terms 
can contribute. Instead, one aims to reduce the errors 
to a small enough level that no extrapolation is needed. 
Second, to estimate the residual errors requires, at the 
least, calculating the full one- loop contributions to im- 
provement coefficients so as to determine their effect on 
the final results. It will be difficult to carry out such cal- 
culations given the complexity of the improved action. 
My major concern is for the calculation of matrix ele- 
ments, for which the perturbative uncertainty enters at 
leading order. For this reason I am skeptical that one will 
be able, in general, to push much beyond a m i n = 0.1 fm 
and obtain precise results for matrix elements. Hopefully 
I will be proven wrong. 

Another approach to improvement is to determine 
a good approximation to the "perfect action" — i.e. an 
action without any discretization errors. This has been 
done at the classical level, but it has proven more diffi- 
cult to find a truncation appropriate for numerical sim- 
ulations which maintains the good properties of the full 
action. For a recent review see Ref.E3. 

4-3 Glueball masses in SU(3) Yang-Mills theory 

The major motivation for improvement is to make simu- 
lations of unquenched QCD tractable. Before discussing 
the status of such simulations, however, I want to give 
two examples from the pure gluon theory for which im- 
provement has allowed substantial progress. 

If one removes not only the sea quarks but also the 
valence quarks from QCD one is left with a pure gauge 
theory with group SU(3) . Unlike quenched QCD, this is 
a physical quantum theory with a unitary S-matrix. It 
is an interesting laboratory for lattice methods because, 
like QCD, its low-energy dynamics are non-perturbative. 
The spectrum consists of interacting glueballs which, if 
massive enough, decay into multi-glueball states. What 
I illustrate here is the progress that has been made in 
studying this theory through the use of improved actions. 

I show in Fig. the latest results, due to Morningstar 
and PeardonE3. They have obtained the masses of the 
lightest 15 glueballs with small errors, including the ex- 
trapolation to the continuum limit. On the LH scale, 
the masses are expressed in terms of a distance scale, 
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Figure 7: Continuum glueball spectrum in SU (3) pure Yang-Mills 
theory. Boxes indicate la errors. Hollow boxes indicate ambiguous 
spin interpretation or ambiguous continuum extrapolation. 



ro, whichjis determined from the quenched heavy quark 
potential EB (and is thus defined in the pure gauge the- 
ory). I stress that these results are definitive numbers 
for dimensionless ratios in the pure gauge theory there 
are no remaining approximations to be dealt with. The 
effects of virtual glueball pairs are fully included. 

To obtain these results two technical tricks were com- 
bined: improvement and anisotropic lattices. The parts 
of the action connecting degrees of freedom in spatial di- 
rections were improved, allowing the use of a large spatial 
lattice spacing: a s = 0.2 — 0.4 fm. The temporal lattice 
spacing was kept smaller, a s /at = 3 — 5, so that the im- 
provement of those parts of the action involving temporal 
separations was not required. This has two advantages. 
First, using a smaller at allows one to trace out the ex- 
ponential fall-off of correlators [C(r) oc exp(— Mr)] more 
accurately and thus extract more precise values for glue- 
ball masses. Second, if one had to improve the action in 
the temporal direction this would introduce unphysical 
"ghost" states which would interfere with the extraction 
of the masses. The use of a larger a s not only reduces the 
computer time required (the calculations have been done 
on DEC- ALPHA workstations), but also makes it easier 
to construct operators having a larger overlap with the 
glueballs. 
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These results are a considerable advance over pre- 
vious work which used the unimproved Wilson gauge 
action E3fEl, and needed to use lattice spacings as small 
as a = 0.05 fm in order to extrapolate to the continuum 
limit. This demonstrates the potential power of improve- 
ment: the minimal lattice spacing is increased by a factor 
of four. In particular, since the improvement was done 
perturbatively, it is a test case for the adequacy of per- 
turbative as opposed to non-perturbative improvement. 
The errors of 0(a 2 ) in the unimproved action are reduced 
to a,g, a s a 2 s and a 2 , and one might be concerned that the 
continuum extrapolations would be difficult if all three 
terms were comparable. In practice, it turns out, with 
one exception, that the a* errors are dominant and con- 
tinuum extrapolations seem straightforward. The excep- 
tion is the scalar glueball (the lightest state), for which 
it appears that the a s a 2 errors are significant, This leads 
to a more difficult continuum extrapolation, and explains 
why the errors for this state in Fig. ^are larger than those 
for most of the more massive states, contrary to the usual 
pattern. I should stress that improvement still works — 
the a s a 2 errors are smaller. than the a 2 errors with unim- 
proved Wilson fermions E3. On the other hand, if one 
wants to reduce the error on the scalar glueball mass 
to a few percent using the improved action one would 
have to work at smaller values of a s than 0.2 fm. This is 
one of the reasons why I chose the conservative estimate 
a mi „ = 0.1 fm for QCD.0 

What are the implications of these results for QCD? 
The masses can be converted to physical units by match- 
ing, say, the quenched rho mass or the heavy-quark po- 
tential to experiment. This is done on the RH scale in 
Fig. ^. Let me stress, however, that the uncertainty in 
these predictions will generically be larger than the 5- 
10% discrepancies between the quenched light meson and 
baryon masses and experiment. This is because when we 
"unquench" glueballs, they have many more decay chan- 
nels than the light mesons and baryons, and, further- 
more, most of the glueballs can mix with quark states. 
Shifts in the masses will, however, be very state depen- 
dent; for example, they will probably be smaller for the 
quark model exotics with J PC = 2 H and H which 
have no quark states to mix with. Because of the po- 
tentially larger quenching errors, I view the results in 
the Figure as only a rough guide to experimental masses, 
and identification with observed glueball candidates dif- 
ficult. I should note that Weingarten and collaborators 
have studied the effects of mixing and decays within 



It is possible to make alterations to the action «oi as to reduce the 
discretization errors in the scalar glueball massca, and this might 
lead to an action for which extrapolations could be done using 
a m i n «s 0.2 fm. While this is an interesting approach, it is not a 
substitute for a systematic improvement program since one must 
make the alterations on a case-by-case basis. 



the quenched approximation, rand conclude that they are 
small for the scalar glueball l3. These conclusions are, 
however, controversialE3. 

Irrespective of its phenomcnological implications, the 
numerically determined glueball spectrum provides a 
benchmark against which to test theoretical models. An 
interesting example is the prediction that states should 
be ordered roughly according to the minimal dimension 
of the operators which create them E3C3. This proposal 
works well for the states in Fig. |7|, which are created by 
operators of dimension 4, 5 and 6. This comparisop_.as 
well as that to the bag model are described by KutiEa. 

4-4 Heavy hybrids 

In pure gauge theory one can introduce an infinitely 
heavy quark and antiquark pair, and calculate the energy 
of the induced gluonic flux as a function of the separation 
r. One can excite the gluonic flux into states classified by 
their angular momentum along the axis and by other dis- 
crete symmetries. There is thus a tower of "potentials" 
which, up to an overall constant, are physical quantities 
in the pure-gauge theory. 

The interest in these potentials is twofold. First, they 
give us insight into the physics of confinement. In the 
standard picture of a flux tube of constant transverse size, 
all the potentials should grow linearly at large r, with the 
excitations corresponding to excited transverse modes. 
At distances of order the transverse size, r ~ 0.5 fm, this 
picture breaks down, and the detailed form of the po- 
tentials contains information about how this happens. 
The second interest is phenomenological. The potentials 
can be used, in the Born-Oppenheimer approximation, to 
calculate the spectrum of quarkonium states. The lowest 
lying potential gives rise to the usual QQ states, while 
the excited potentials correspond to hybrid QQG states. 
A prediction for the mass of the lightest hybrids in say, 
the bb system would be very interesting. Of course, such 
a prediction would be in the quenched approximation, 
since the effects of virtual quarks (light and heavy) have 
been ignored. 

Recent work by Juge, Kuti and Morningstar has coia- 
siderably extended our knowledge of these potentials Ell. 
These authors use the same improved, anisotropic lat- 
tice action as in the glueball calculation of Ref. E3, with 
spatial and temporal lattice spacings in the ranges a s — 
0.12 - 0.29 fm and a t = 0.04 - 0.06 fm respectively. Com- 
pared to previous work using the unimproved action E3, 
they are able to calculate more excited potentials and to 
extend the results out to much greater separations. A 
subset of the new results are shown in Fig. @. For clarity, 
the potentials are shown only out to about 2 fm, although 
they have been determined out to about 4 fm. The group 
theoretic symbols indicate different representations of the 
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r/(2r ) 

Figure 8: Potentials in iqwenched QCD after continuum extrapola- 
tion (data of Juge et a/.EJ plotted in Ref.Ej). Units are ro ~ 0.5 fm. 
The two lines for each curve indicate the estimated error. In the 
effective string theory the first four predicted bands are: (I)£ 9 ; (II) 
II U ; (III) £' + , n 9 and A 9 ; (IV) E±, II ^ and A„. 

gluon flux, the details of which can be found in the orig- 
inal work. 

What I want to stress is that, as for the glueball 
spectrum, numerical methods have produced detailed, 
accurate results for non-perturbative quantities. These 
results are good enough to test theoretical models. For 
example, at long distances, an effective string theory de- 
scription should apply, which predicts that the excited 
states should fall into bands. This band structure is in- 
deed observed by about 2 fm and is visible in the fig- 
ure. The string model also predicts that the differences 
between the levels asymptotically approach multiples of 
7r/r, but this turns out not to workEll. This discrep- 
ancy may be due to a flaw in the model itself, as dis- 
cussed by.Ref.t3, or to underestimated systematic errors 
(see Ref. lB for a discussion pf possibilities). An impor- 
tant point, stressed by KutiEa, is that the fact that the 
ground state potential remains very close to linear all 
the way down to r ss 0.5 fm does not necessarily provide 
evidence in favor of the sting picture at these distances. 
Such "precocious" linearity is predicted by models (e.g. 
the adiabatic bag modelE3) to occur well before the onset 
of string-like behavior. 

At intermediate distances, r w 1 fm, the excited 



states undergo a complicated pattern of level crossings. 
While some of this behavior can be understood in teems 
of degeneracies that are predicted for r — (Ref. Ej), 
it remains a challenge to understand these levels in de- 
tail. This work of interpretation is just beginning — for 
the moment the important progress is having data good 
enough to compare with models. 

Let me turn now to the predictions for the masses of 
hybrids. In the Born-Oppenheimer aunroximation, the 
lightest bbg level is eightfold degenerateE2l. This level con- 
tains a J = 1 , which could be produced in e + e~ 
annihilation, and three exotics with J PC = 1 h , H , 
and 2" 1 . This degeneracy is lifted when one includes 
retardation effects, as can be done using non-relativistic 
heavy quarks simulated on the lattice. Such calculations 
have been done by Manke et a/. lid and by Juge et al.E3. 
The retardation effects are found to be small — for exam- 
ple, all eight levels remain close in mass. The detailed 
predictions from the two groups are in slight disagree- 
ment. For example, the mass of the bbg hybrid_with 
J PC = 1-+ is fooind to be 11.1(1) GeV by Ref.O and 
10.8 GeV by Ref.Ej, both in the quenched approximation. 
This is an interesting mass range for hybrids, because, al- 
though above the BB threshold at 10.56 GeV, it brackets 
the S-wave+P-wave B}B threshold at 11.01 GeV. There 
are theoretical arguments suggesting that the dominant 
decay mechanism of a hybrid is to an S- and a P-wave 
meson i3. Thus, if the hybrid lies below or close to the 
BjB threshold, it could be a narrow state, and thus ob- 
servable. Improvements in the quenched results are likely 
in the next year. 

Calculations of ccg irybrids have been done by 
the MILC collaboration La, using improved relativis- 
tic quacks, and by Juge et al., using non-relativistic 
quarksN. MILC finds 4.3+^ GeV for the 1~+ state, 
and the result from Juge et al is consistent with this. 

5 Simulations of two-flavor "QCD" 

An important step on the way to simulating full QCD is 
the study of the theory with two flavors of light quarks. 
As noted in sec. |[ the use of an improved fermion action 
will likely be essential in order to simulate this theory 
with quarks light enough to allow a controlled extrapo- 
lation to the chiral limit. The hope is that improvement 
will work well enough to allow a continuum extrapolation 
to be done using lattice spacings no smaller than 0.1 fm. 
The successes of improvement in the pure-gauge theory 
auger well for this program. 

Machines sustaining several hundred Gigaflops have 
been used to do preliminary studies of two-flavor QCD. 
The extra cost of including dynamical fermions restricts 
these simulations to larger lattice spacings (a > 0.1 fm), 
smaller volumes (L < 2.4 fm), and larger quark masses 
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Figure 9: CP-PACS results for mx* and versus lattice spacing 
in two-flavor QCD (filled symbols), compared to quenched results 
(open symbols including continuum extrapolation). The strange 
quark mass is fixed by setting mx to its experimental value. All 
results are shown after chiral extrapolation, which is done using 
a polynomial fit for the two-flavor data, and a quenched chiral 
perturbation theory fit for the quenched data. The diamonds at 
a = show the experimental values. 



(m^/mp > 0.6) than the state-of-the-art quenched simu- 
lations discussed above. The most important of these re- 
strictions is on the quark masses — reliable extrapolations 
to physical quark masses are not yet possible. Neverthe- 
less, some interesting results have been obtained and I 
will highlight a few of these. 

Several collaborations have undertaken simulations 
aimed at studying the spectaum of light hadrons in 
two-flavor QCD: CP-PACS eMB, using a perturbatively 
0(a) improved fermion action and an improved gauge 
action-determined by renormalization group methods; 
MILCL3, using staggered fermkms (which are automati- 
cally 0(a) improved); SESAMO, using unimproved Wil- 
son fermions; and UKQCD L3, using non-perturbatively 
0(a) improved Wilson fermions. Of these, the CP-PACS 
and MILC studies are the most extensive, using several 
lattice spacings and the widest ranges of quark masses. 

An interesting first question to study is whether the 
previously noted discrepancy between the hyperfmc split- 
ting in quenched QCD and experiment is reduced in two- 
flavor QCD. The preliminary CP-PACS results, shown 
in Fig. |[ do show a significant increase in this splitting 
which moves the results toward the experimental value. 
This is encouraging, although much more work is needed 
to control the chiral and continuum extrapolations. 

Other quantities for which a significant effect of un- 
quenching has been observed are the quark masses (dis- 



cussed in the next section), and the heavy quark poten- 
tial. The latter should behave differently both at short 
and at long distances. At short distances, the one-gluon 
exchange contribution oc a s (r) / r should be stronger in 
the unquenched theory because the coupling runs more 
slowly due to the dynamical quarks. This effect has been 
observed with the expected magnitudes. 

At long distances a more striking phenomenon 
should occur — the string should break! This should oc- 
cur at r ~ 1.5 fm, when it becomes energetically favor- 
able for the flux tube to decay into two sources screened 
by light (anti)quarks. At this point the potential should 
bend over and reach a plateau. Surprisingly, this rather 
basic phenomenon has not yet been seen in simulations. 
A year ago, as summarized in a review by Giisken l3, 
results for the potential had been extended out to dis- 
tances at which the expected flattening should begin. 
This year, CP-PACS has pushed the calculation of the 
potential out to somewhat longer distances, yet still finds 
that the potential keeps on rising La. There has been 
progress, however., in understanding why string-breaking 
is not observedE^H. The key point is that the operators 
used to create the gluon flux have very small ove 
with the broken-string state. The authors of Refs. EiH 
find that if they include in their bases operators designed 
to have substantial overlap with the broken-string state, 
then they clearly observe string breaking. These studies 
are done, however, on theories with scalar quarks, and it 
would be very interesting to repeat them on QCD itself.H 



6 Hadronic matrix elements 

A major reason for simulating lattice QCD is to extract 
results for the hadronic matrix elements of various opera- 
tors. We would like to have both the numerical values of 
these matrix elements and reliable estimates of the uncer- 
tainties in these values. In this section I report on results 
for those quantities for which the calculations are mature 
enough to include a continuum extrapolation within the 
quenched approximation. These results are collected in 
Table || and I will discuss them in more detail in the fol- 
lowing. My aim is to allow the non-specialist to make a 
judgement on the reliability of the error estimates. For 
more technical details,.-!, refer the reader to the reviews 
of MartinelliEJ, GuptaEJ and LubiczEJ (mainly for light 
quark matrix-elements), and of Flynn and SachrajdaEZI 
and Draper E3 (for heavy- light matrix elements). These 
reviews also contain updates on many other interesting 
quantities for which the calculations are less developed. 



s Stri»g| breaking has already been observed in finite temperature 
QCDEj. Here one has the option of using the correlation between 
Polyakov lines to measure the potential and these have a good 
overlap with the broken-string state as well. 
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Table 3: Summary of results for matrix elements. The third column contains estimates of the effect of removing the quenched approximation. 
Quark masses are quoted in the MS scheme at a scale of 2GeV. All B— parameters quoted are renormalization group invariant, defined 
including non-leading order terms. Results without references are my estimates. 



Quantity 


Quenched Result 


Effect of 
unquenching 


Comments 


m = 0.5(m u +md) 


4-5 MeV 


« 40% reduction 


Non-perturbative matching crucial 


m s 


100 - 150 MeV 


« 40% reduction EH 


Using ttik and m§ gives different results 


B K 


0.85 ± 0.06 (Ref. 


103 \ 


x 1.10 ±0.16 (Ref.@) 


Unquenched result is 0.94 ± 0.06 ± 0.14 


Ib 


165 ±20 MeV (Rcf.N) 


+2 n 3 MeV (Ref.N) 


Ref. 104 quotes 175 ± 25 MeV (quenched) 


B b 


1.4 ±0.1 (Ref.l 




See text 




Ib \J B b 


195 ±25 MeV 


See text 


Ref.El quotes 201(42) MeV (quenched) 


' V &Bb 


1.14 1* (Ref.il) 


See text 


Most reliable of all B-meson results 


Id, 


220+^ MeV (Ref.N) 


t\ 7 MeV (Ref.N) 


Compare with experiment 



Three other quantities, a s , m c and mb, are also well 
determined from lattice simulations. I do not discuss 
them here because there has been little new work on 
them in the last year. Detailed discussions are ghzrn in 
earlier reviews: for a s see the particle data tablesLJ-and 
references therein; for heavy quark masses see Ref.lHa. 

6.1 Light quark masses 

Quark masses are fundamental parameters in the stan- 
dard model. In principle, the masses of the up, down 
and strange quarks are encoded in the masses of the 
light hadrons, but to determine them requires a non- 
perturbative method. Lattice calculations will ultimately 
provide results with all errors under control, but at 
present they are plagued by the uncertainties due to the 
quenched or partially quenched approximations. Calcu- 
lations using QCD sum rules remain for the moment of 
comparable reliability, and competition between the two 
"camps" has produced useful reanalyses on both sides. 

Quark masses are also interesting for a more practi- 
cal reason. In the usual parameterization of the standard 
model prediction for e'/e in terms of B— parameters, the 
dominant contributions are inversely proportional to m 2 s . 
If, as lattice calculations suggest, m s is smaller than ex- 
pected a few years ago, then the prediction for e'/e is 
likely to increase. Unfortunately, this is not a generic 
statement because there_axe delicate cancelations. For 
recent reviews see Refs.E3o. 

The extraction of masses from lattice calculations 



is, in principle, straightforward. Ignoring isospin break- 
ing, one has two dials to turn: the bare average light 
quark mass m and the bare strange quark mass m s . One 
adjusts these until the results for two mass ratios, say 
m^/irip and mx/m^, agree with experiment. One then 
determines the "physical" quark mass (defined in some 
renormalization scheme at a fixed scale) by matching the 
lattice and continuum definitions of quark masses. There 
are several ways of doing this — but the key point here is 
that all require knowledge of matching factors relating 
lattice and continuum operators. For example, one can 
use the ratio of the divergence of the axial current to the 
pseudoscalar density to define what is called the "axial 
Ward identity" mass: 

2 _ cont = (o|^ ont l g ) 

(0|F cont |i?) 

z A (Q[ yijf Iff) 

Zp (0|P lat |#) ' [ ' 

Here \H) is an arbitrary hadronic state (usually the 
pion), "cont" and "lat" refer respectively to continuum 
and lattice regularized operators, and the matrix ele- 
ments are to be evaluated with the bare quark masses 
set to their previously determined values. In this method 
the matching factors one needs are Za and Zp. 

Let me first discuss the average of the up and down 
quark masses, m. The state-of-the-art a year ago was, for 
the most part, to use one-loop perturbative results for 
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Table 4: Comparison of perturbative and non-perturbative match- 
ing factors. For WBson fermions the pseudoscalar density is defined 
in the RI scheme E3, while for staggered fermions the MS scheme 
is used. Both are defined at a scale fi = 2 GeV. All results are for 
quenched QCD with /3 = 6/g 2 = 6 (a ss 0.1 fm). 



Fermion type 


Quantity 


One-loop 


Non-pert. 


Wilson y 


Z A 


0.78 


0.81(1) 


Z P 


0.62 


0.45(1) 


Non-pert, 
improved. 
Wilson O 


z P 


0.59 


0.39(3) 


StaggeredS 


Z A 


1.0 


1.0 


z P 


0.54 


0.34 



matching factors, and to do chiral extrapolations with 
the lowest order polynomial that gave a good fit. Af- 
ter continuum extrapolation of all available data Bhat- 
tacharya and Gupta concluded that, in the quenched 
approximation, 

m = 4.1 ±0.1 MeV (Wilson) 

= 3.8 ±0.1 MeV (Perturbatively Improved Wilson) 
= 3.5 ± 0.1 MeV (staggered) 

for the different fermion discretizations. I emphasize that 
these results are for an MS scale of fj, = 2 GeV (and 
not 1 GeV as often used in the sum-rule literature). Re- 
sults using different discretizations of the fermion action 
should agree in the continuum limit. The discrepancies 
between them were attributed to systematic errors asso- 
ciated with uncertainties in the continuum extrapolations 
and matching factors. Other groups, making different 
assumptions about these systematics, found somewhat 
larger values for the quark masses. 

The most important improvement in the last year 
is that matching factors have been determined non- 
perturbatively. The technology for doing this has been 
developed over the last decade, and is an interesting ex- 
ample of "applied field theory" . The normalization of 
the axial current is determined by enforcing Ward iden- 
tities E3, while that of the pseudoscalar density (and 
other divergent operators) is_d£jtermined using off-shell 
quark correlation functions ol3. I give several exam- 
ples of the difference between one-loop perturbative and 
non-perturbative results in Table ^ The one-loop result 
works well for Za but not for Zp, which is not unex- 
pected given the sizes of the one-loop corrections in the 
two cases. This is an example of the importance of us- 
ing non-perturbative methods to normalize operators — a 



point I stressed previously when discussing the improve- 
ment program. Another success of non-perturbative nor- 
malization is that it brings results for quark masses deter- 
mined using different methods (vector versus axial Ward 
Identity) into agreement E£j. 

Referring back to Eq. O, we see that smaller val- 
ues for Zp imply larger quark masses. For example, the 
new result for staggered fermions, m = 4.2(3) MeV lB, 
is 20% larger than that quoted above. The inclusion 
of non-perturbative matching factors does not, how- 
ever, seem to bring the results using different types of 
fermion into closer agreement. For example, there is, for 
the first time, a calculation of quark masses using non- 
perturbatively O(a) improved Wilson fermions ES, which 
gives m = 4.9(4) MeV. Although this result does not in- 
clude a continuum extrapolation it is from a small lattice 
spacing, a w 0.07 fm, for which the 0(a 2 ) discretization 
errors are likely to be small. 

Another advance in the last year concerns the chi- 
ral fits. As discussed above, there is reasonable evidence 
in the quenched spectrum for the artifacts predicted by 
quenched chiral perturbation theory. Fitting including 
these artifacts results in quark masses that, at finite lat- 
tice spacing, are significantly reduced compared to those 
from polynomial fitsE3-Eil. It turns out, however, that, 
with present statistics, the shift is not significant after 
continuum extrapolation. 

All in all, despite the progress that has been made, 
we are lacking the comprehensive data sets, including 
non-perturbative matching factors, which would allow 
fully controlled chiral and continuum extrapolations. My 
summary of the present situation is that, in the quenched 
approximation, m lies in the range 4-5 MeV.0 I expect 
some clarification in the upcoming year, particularly us- 
ing the matching factors calculated using the Schrodinger 
functional E3. 

The upwards movement of the quenched results 
brings them into reasonably comfortable consistent: 
with the bounds derived using QCD sum rules | 101 H 1CL 
But what are the quenching errors? Early calculations 
suggested that quark masses were smaller when quark 
loops were included, and_this has been confirmed by the 
results from CP-PACS LJ. Their preliminary analysis 
suggests that two flavors of dynamical quarks reduce the 
light quark masses by about 40%, which might well bring 
them into contradiction with the sum rule bounds. It will 
be very interesting to see how things evolve over the next 
few years. There must be unexpectedly large systematic 
errors in cither the lattice or the sum rule methods (or 

both) Suggestions for the latter have been discussed in 

Ref. 
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h In the talk I mistakenly quoted a larger range 4 — 6 MeV. The 
upper end of this range was based on the result from Ref. £3 in the 
RI, rather than the MS, scheme. 
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The systematic errors in determining m occur also in 
the determination of the strange quark mass. In addition, 
as is clear from comparing Figs. ^ and||, a different value 
results if one uses the tuk or to fix the strange mass]] 
Using rriK leads to m s /m ~ 25 if the non-linearities in 
m^/rriq are small, whereas using turns out to give 
m s /m w 30. Thus I have converted m = 4 — 5MeV to 
m s = 100 — 150 MeV for my estimate of the likely range 
for the quenched strange quark mass. 

6.2 The kaon B— parameter 

Little has changed in the_result for Bk since the reviews 
by Flynn u and myself E3. Nevertheless, since it is an 
important number for phenomenology, I want to stress a 
few points. 

The most reliable determination comes from stag- 
gered fepnions, and the best data by far is that from 



JLQCD 



103 



In a tour de force they have distinguished 



between discretization errors and errors due to the use 
of one-loop perturbative expressions for matching fac- 
tors. The former are proportional to a 2 while the latter 
vary as [a s (a)] 2 oc [In a] -2 . Their final result is that 
-Bfs-(NDR, /i = 2GeV) = 0.628(42). The relatively large 
error is almost entirely due to the uncertainties in fit- 
ting the a dependence their statistical errors are much 
smaller. This is yet another indication that accurate 
numbers for matrix elements require non-perturbative 
determinations of matching factors. 

The JLQCD data also contains another warning. Al- 
though the calculation is automatically improved, with 
the discretization errors being of 0(a 2 ), the contin- 
uum extrapolation requires accurate data down to a = 
0.05 fm. This is because the dimensionful coefficients 
multiplying powers of a turn out to be large. In the 
notation of Eq. B even though Ai = 0, A2 is a few times 
larger than Aqcd, as are the higher order A„. Thus im- 
provement does not always work as expected. 

Converting the JLQCD result to the renormalization 
group, invariant quantity Bk (using the assumptions of 
Ref.EJ) gives the result quoted in Table ||. Fortunately, as 
noted by GuptaEJ, the result for Bk is quite insensitive 
to whether this conversion is done using the quenched or 
unquenched formula. This final quenched result differs 
from that quoted in 1996 because the final JLQCD num- 
ber is slightly higher than their preliminary result due to 
the more sophisticated continuum extrapolation. 

The estimate of quenching errors has not changed 
significantly since 1996 Numerical data from un- 



1 There is an analogous dependence of m on whether it is determined 
using m w /m p or m w /m]y, but, as can be seen from Fig. |i| the 
difference is less statistically significant than that for m a . It is 
conventional to fix m using m p rather than m jv , because the former 
has a smaller statistical error. 



quenched simulations with two and four flavors implies 
that quark loops increase Bk by a factor of 1.05 ± 
0.15, where the error estimate is quite conservative. 
The change from degenerate quarks to physical non- 
degenerate quarks in the kaon is estimated, using chiral 
perturbation theory, to increase Bk by another factor of 
1.05 ± 0.05. Here the error estimate is less conservative. 
Combining these two factors gives 1.10 ±0.16, where the 
errors have been combined in quadrature. Multiplying 
the quenched result by this factor gives the final result 
quoted in the table. Whether one should prefer this num- 
ber to the quenched result, supplemented by an estimate 
of the quenching error, is a matter of taste. I prefer to 
include the estimates of the effect of unquenching, along 
with their fairly conservative errors, since they are based 
on reasonably reliable arguments. 

6.3 B— meson matrix elements 

The study of heavy-light mesons and baryons is now the 
major focus of that part of the lattice community ori- 
ented towards phenomenology. This is certainly appro- 
priate given that B-factories promise a wealth of pre- 
cise data in the near future. Lattice calculations will 
likely play an important role in the extraction of the un- 
derlying parameters of the standard model from these 
data. Although I report here results only for decay con- 
stants and B— parameters, I want to stress that fully con- 
trolled quenched results for many other quantities are 
likely to appear soon. Examples include form factors for 
B -> D^lv, B -> {p,ir)lv, D -> K^tv and B -> AT* 7 , 
and the spectrum of excited hadrons containing b quarks. 
For further details see the above-mentioned reviewsEZIH. 

The most important theoretical issue in heavy quark 
calculations on the lattice is the accuracy with which 
one can simulate the b quark. Direct simulation of a 
fo-quark with a relativistic lattice action is not possible 
because discretization errors are large for present lat- 
tice spacings since m^a 3> 1. To avoid these errors, 
one can either use a relativistic action but extrapolate 
from lighter quarks (typically somewhat heavier than 
the charm quark), or treat the 6— quark using a non- 
relativistic expansion. Both approaches have potentially 
large systematic errors. In the relativistic approach these 
are from discretization errors, which can be reduced using 
improved actions and operators. In the non-relativistic 
approach the errors arise from the truncation of the non- 
relativistic expansion, and the need (so far at least) to 
match the non-relativistic theory with QCD using per- 
turbation theory truncated to one-loop order. Although 
the various groups are making quite sophisticated error 
estimates for their respective calculations, the uncertain- 
ties in the errors themselves are large. | This makes it 

•>For a clear discussion of the errors in the relativistic and non- 
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UKQCD98a 
UKQCD98b 
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Figure 10: Comparison of results for Jb, taken from Draper E3. 
Statistical and systematic errors are superimposed, the former in- 
dicated by the larger tick marks. The quoted average is also shown, 
together with the results included in determining it. 



difficult, at least for someone not involved in the calcula- 
tions, to make a meaningful average of results from the 
different groups, and so I haiffi chosen to quote averages 
from the review of Draper E3, iL available, or from the 
review of Flynn and SachrajdaEH. Fortunately, in most 
cases, the results from different methods are consistent, 
so the central values are not likely to change much with 
a different weighting of the different results. I have il- 
lustrated this robustness by quoting examples of results 
from other reviews in the table. 

For the decays constants, Draper averages results 
from the following collaborations: GLOK t23, who jasc 
non-relativistic heavy quarks; FNALtL°a and JLQCDt23, 
who use improved Wilson —quarks interpreted non- 
relativistically following Ref. MILCcj, who extrap- 
olate using Jioimproved W-ilion quarks, interpreted fol- 
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who use both unimproved 
The consistency between 



lowing Ref. and APE 
and improved Wilson quarks 
these results is illustrated in Fig. [l0[ along with Draper's 
average. Also shown are various numbers which are not 
included in the average: old results which do not have 
as good control over systematic errors, and new results 



from various collaborations, which contain various im- 
provements, but which are preliminary. 

The quantities of particular interest to phenomenol- 
ogy are the last three in Table |^: f B \J~Bb~ is needed to 
extract information on the CKM elements from B — B 
mixing; £ is needed to determine the constraints from 
B s — B s mixing; and fr> s can be compared with exper- 
imental results. The errors in the quenched results for 
these quantities are relatively small, and it is now crucial 
to obtain estimates of the quenching error. For the de- 
cay constants themselves, such estimates have been made 
using numerical simulations by the MILG-collaboration, 
and are included in the table. Reference Ii2a gives a de- 
tailed exposition of how these estimates are made, and 
notes that they are quite uncertain — the error on the 
quenching error is itself estimated to be ~ 50%. For 
the B parameters there is little numerical guidance as 
to the size of quenching errors. One can make theo- 
retical estimates, based ,uaaLL comparing chiral loops in 



These suggest quenching 



LIS , 33 



relativistic approaches see Ref.li£j and Ref.li£j, respectively 



quenched and full QCD 
errors as large as ~ 10% in the individual B parameters, 
though somewhat smaller in the ratio BbJBb- Given 
the roughness of these estimates, as well as the uncer- 
tainty in the quenching errors for the decay constants, 
it is probably premature to attempt an estimate of the 

unquenched results for Jb\J Bb and £. Nevertheless, one 
requires numbers in order to place bounds on CKM pa- 
rameters. My best estimates for the unquenched num- 
bers are 1b\[b~ b = 210 ± 40 MeV and f = 1.14(13). 
To obtain these I have included the MILC estimate for 
the shift in fs due to unquenching, but inflated their 
quenching errors by 50%, and taken quenching errors on 
B B , IbJ/b, and B B JB B to be 10%, 10% and 5% re- 
spectively, based on the chiral loop estimates. Draper 
has given similar unquenched estimatesE3. 

7 Exactly massless quarks on the lattice 

My final topic concerns discretizations of the fermion ac- 
tion which exactly preserve (or almost preserve) the full 
global chiral symmetries of continuum QCD at non-zero 
lattice spacing. This is to be contrasted with existing dis- 
cretizations (Wilson and staggered fermions), for which 
one cannot separate the chiral and continuum limits — 
chiral symmetry is only fully restored in the continuum 
limit. Separating these limits is very attractive from both 
a formal and a practical point of view. Continuum discus- 
sions of chiral symmetry breaking, topological properties, 
and rigorous inequalities often rely on exact chiral sym- 
metry in the massless limit. These arguments carry over 
directly onto the lattice if chiral symmetry is maintained. 
Practical calculations of quantities constrained by chiral 
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symmetry, such as the quark condensate and kaon mixing 
and decay amplitudes would also be greatly simplified if 
there were no explicit chiral symmetry breaking in the 
lattice action. During the last year there has been a 
flurry of activity on this possibility, and I will attempt to 
sketch the essence of this work and its implications. For 
more details see the recent review of NiedermayerS, and 
the contribution by Neubergerla. 

Recall that, formally, chiral symmetry is restored in 
the massless limit because the Dirac operator anticom- 
mutes with 75: {^,75} = 0. When one discretizes the 
fermion action, and removes the doublers in a reasonable 
way (local Dirac operator with correct continuum limit), 
one finds that the anticommutator no longer vanishes 



{Plat, 75} =0(6 



(12) 



This result follows from the famous Nielsen-Ninomiya no- 
go theorem 1113. Two lines of work have attempted to 
evade the explicit chiral symmetry breaking on the RHS 
of Eq. [L2UThe first was begun by Ginsparg and Wilson 
in 1982E1 They noted that the 0(a) term in Eq. [TJ 
could, in principle, be made to vanish in physical matrix 
elements by a suitable choice of fermion action. They 
formulated a condition, now known as the "Ginsparg- 
Wilson relation" , the simplest form of which is 



{Plat,75} = aBlat75p 



lat 



(13) 



If this relation is satisfied, then chiral Ward identities 
are violated only by contact terms (because the two fac- 
tors of D on the rhs of Eq. [ll] cancel the propagators 
joining the lattice artifact to other parts of correlation 
functions), and thus all physical consequences of these 
identities remain valid.^ In other words, the Ginsparg- 
Wilson relation ensures that there is, in effect, a chiral 
symmetry. Despite these nice properties, no explicit solu- 
tions to the relation were found when gauge interactions 
are included. The difficultly is clear from an equivalent 
form of the relation, {^^,75} = a 75, which requires 
the anticommutator of the inverse Dirac operator with 
75 to be local. Thus the Ginsparg- Wilson relation lay in 
obscurity for many years. 

The second line of work was the attempt, by many 
authors over many years, to come up with a discretized 
version of a chiral gauge theory (or at least a lattice the- 
ory which one could show became a chiral gauge theory 
in the continuum limit). A very interesting approach 
to this problem was initiated by Kaplan's domain wall 
proposal EM As explained in Ref. □, this and other 



work prompted Narayanan and Neuberger to develop 
the "overlap" formulation of chiral gauge theories E 17 



k Ginsparg and Wilson did not come upon their relation in this way, 
but rather used a renormalization group strgiiment. This interpre- 
tation of their result is due to Hasenfratzlii^ . 



Whether this proposal works for chiral gauge theories 
is a matter of some controversy, but it is uncontroversial 
that the "overlap" provides a formulation of vector gauge 
theories such as QCD with exactly massless fermions. It 
evades the no-go theorem by having, in effect, an infinite 
number of regulator fields. This solution to the prob- 
lem of exactly massless quarks in vector theories also 
lay in relative obscurity for the last few years, perhaps 
because a practical implementation of the overlap formu- 
lation seemed difficult. 

In the last two years these threads have come to- 
gether. Neuberger showed how one could integrate out 
the extra degrees of freedom in the overlap formulation 
and obtain a relatively simple form for a lattice Dirac op- 
erator corresponding to exactly massless quarksEM This 
he called the "overlap-Dirac operator" . The second step 
was the "rediscovery" of the Ginsparg- Wilson relation: 



^classically perfect" fermion 
1 1 These actions (reviewed 



Peter Hascnfratz noted tha 
actions-satisfy the relation 
in Ref.E3) thus provide another solution to the problem of 
discretizing a vector gauge theory with massless quarks. 
Neuberger then noted that the overlap-Dicac operator 
also satisfied the Ginsparg- Wilson relationllHS. 

Finally, Liischcr clarified the previous work by 
showing that any lattice Dirac operator satisfying the 
Ginsparg- Wilson_cclation has an exact, albeit non-local, 



chiral symmetryEziJ. Furthermore, he showed with a very 
compact argument how "Ginsparg- Wilson fermions" sat- 
isfy an exact lattice version of the continuum Atiyah- 
Singer index theorem. His work has led to further devel- 



122 



opments in the formulation of chiral gauge theories^ 

The end result is that we now have explicit forms of 
lattice fermion actions with exact chiral symmetry, and 
a simple and general formalism for understanding them 
based on the Ginsparg- Wilson relation. This advance has 
spawned a flood both of theoretical work aimed at deter- 
mining the consequences of the Ginsparg- Wilson relation 
in detail, and at finding a practical implementation. 

7.1 Numerical implementations 

My previous sketch left out what may be the most impor- 
tant development from a practical point of view, namely 
"domain-wall fermions" . Kaplan's proposal for a chiral 
gauge theory requires an extra, fifth dimension of infi- 
nite extent, in which the chiral theory lives on a four- 
dimensional domain-wall. Truncating to a finite fifth di- 
mension introduces a second domain-wall inhabited by 
fermions of opposite chirality, so that overall pne has a 



vector theory. Thus, as stressed by ShamirE£3, domain- 
walls provide a practical method for studying QCD. The 
coupling between the left and right-handed fermions is 
suppressed exponentially by the distance between the 
two walls, and so violations of chiral symmetry are con- 



18 



trollable and can, in principle, be made as small as one 
wants. One has again succeeded in separating the chiral 
and continuum limits. 

In the last year there have been extensive numerical 
simulations using domain- wall fermions, and they indi- 
cate that this is a very promising approach. In particu- 
lar, it may be possible to get away with only 10 — 20 sites 
in the extra dimension. Of course, the simulations slow 
down by this factor, but this may be more than compen- 
sated for by the lack of chiral symmetry breaking, and 
by the fact that the domain-wall action is automatically 
0(a) improved. For more details see the recent review 
There have also been investigations of the 



by Blum El 



properties of the overlap-Dirac opeator and the feasibil- 
ity of numerically inverting it. The non- locality makes 
inversion slow, but the potential gains are great. I ex- 
pect a lot of development in these directions during the 
upcoming year. 
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